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1. Introduction 



The discrete self-trapping (DST) equation was introduced by Eilbeck, Scott and 
Lomdhal fl]] to model the nonlinear dynamics of small molecules, such as ammonia, 
acetylene, benzene, as well as large molecules, such as acetanilide. In simple terms, it 
consists of a set of n nondissipative anharmonic oscillators coupled through dispersive 
interactions. Due to the nonlinearity this system can have complicated dynamical 
behaviour going from quasiperiodic motion to chaos 0, ||. The DST equation is also 
found in connection with physical problems in different areas such as stabilization of 
high frequency vibrations in biomolecular dynamics [Q, arrays of coupled nonlinear 
waveguides in nonlinear optics [^, quasiparticle motion on a dimer 0. In the case of 
two degrees of freedom n = 2 (DST dimer) the system is integrable having, besides 
the Hamiltonian (energy), another conserved quantity, the norm (number of particles 
in the quantum case). The integrability properties of the classical and quantum DST 
dimer were studied in detail by several methods such as the number state method 
[0], the algebraic Bethe ansatz and the method of separation of variables p. For 
more than two degrees of freedoms an integrable case of the DST system was found 
and studied in [10|. This integrable case is close to the Toda lattice and coincides for 



n = 2 with the usual DST dimer. 

The quantum Hamiltonian H of the integrable DST model contains (n+1) parame- 
ters Ci, . . . , c„, 6 and is defined as a second-order differential operator (here di = d/ dxi) 

" /I 1 \ 

^ = E [-j^^l^^i + + -)XA + (1.1) 

acting in the space C[a;] of polynomials of n variables {xi, . . = x. In ( |1 . 1|) and 
other similar formulas we always assume the periodic boundary conditions Xn+i = xi. 
The Hamiltonian H obviously commutes with the number-of-particles operator 

n 

N = J2xA. (1.2) 



i=l 



As shown in section , H and can be included in a commutative ring of differential 
operators generated by a basis of n operators, the fact allowing to claim the quantum 
integrability of the system. 

The multiplication operators Xj and the respective differentiations di can be consid- 
ered as generators of a Heisenberg algebra (creation/annihilation operators). There 
exists a well-known scalar product on C[x] (holomorphic representation) such that 
Xi and di become mutually adjoint dj = Xi. The corresponding Hamiltonian H is 
self-adjoint, however, only in the dimer case n = 2. In the general n > 1 case, no 
involution rendering H self-adjoint is known. The Hilbert space structure is, however, 
quite irrelevant for the kind of problems we are interested in and will be completely 
ignored throughout the paper. 

The DST-chain can be considered as a degenerate case of the Heisenberg magnetic 
chain, though not as degenerate as Toda lattice. It makes the DST-chain a good tool 
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for studying various techniques applicable to integrable models since it requires more 
effort than Toda lattice but still is simpler than the generic magnetic chain. 

The main purpose of the present paper is to construct an analog of Baxter's Q- 



operator |]TT| for the integrable DST model. The Q-operator Qx by definition shares 
the set of eigenvectors with the Hamiltonians Hi, and its eigenvalues are polynomials 
in A satisfying a finite-difference equation known as Baxter or separation equation. 
As was shown in [|12| on the example of the periodic Toda lattice, the similarity 



transformation O QxOQ^ turns in the classical limit into a classical Backlund 
transformation that is a one-parametric family of canonical transformations preserving 
the commuting Hamiltonians. Later, in [13| for the classical Backlund transfomations 
the property of spectrality was described which is the classical counterpart of the 
separation equation for the eigenvalues of Q\. In the present paper we follow the 
approach of [|I^ studying first the classical case and paying the special attention to 
the spectrality property of the corresponding Backlund transformation. 
Our main result (see sections 4-7) is the following integral equation 



/ d^i... / de„ 



1 

n^r(A + l-c,)e-^^(-e, 

j = l 



,Ci~X—l 



i'{---,ykCk + byk+i,---) (1.3) 

(1.4) 



= q{X)ilj{xi,...,Xn), g(A) G C[A], 
for the polynomial eigenf unctions ip E C[x] of the Hamiltonian (|1 . 1|) 

Htjj{xi, ...,Xn) = h'ip{xi, . . . , x„). (1.5) 

The structure of the paper is the following. In section 2 we consider the classical 
version of the integrable DST-chain and describe its relation to the Toda lattice and 
the isotropic Heisenberg magnetic chain. Our construction of Backlund transformati- 
on generalizes well known results for the Toda lattice. Following we also study 
the dual Lax matrix and the corresponding dual Backlund transformation in section 
3. 

In section 4 we discuss the quantization of the integrable DST model and present 
a list of properties of Baxter's Q-operator. In section 5, following the approach of 
[0, we construct a Q-operator Qx for the quantum DST chain as the trace of a 



monodromy matrix with an infinite-dimensional auxiliary space. In the spirit of |jT2 
we consider Qx as an integral operator in C[x] and find in section 6 its kernel and 
contour of integration. In the same section we study analiticity properties of Qx in the 
parameter A, prove that its matrix elements in the monomial basis are polynomials 
in A and give explicit formulas for its action on polynomials. We consider in details 
the simplest n = 1 case where the Q-operator provides an integral representation for 
classical orthogonal polynomials (Charlier polynomials). In section 7 we prove that 
Qx satisfies a finite-difference equation in the parameter A. In the last, 8th section we 
discuss possible generalizations and applications of our results. 
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2. Classical case 



In this section we consider the classical integrable DST-chain ||10[. The model is 
described in terms of n pairs of canonical variables (Xj, Xj), i = 1, . . . , n 

{Xi, Xj} = {xi, Xj} = 0, {Xi, Xj} = Sij (2.1) 

(the periodicity convention Xi+n = Xi, Xi+n = is always assumed for the indices of 
Xi and Xj). 

The canonical momenta Xj replace in the classical case the differential operators (9j. 
As mentioned before, in the quantum case we do not make any assumptions about self- 
adjointness of the observables. Respectively, we allow the classical variables (Xj,Xj) 
to be complex. 

To construct n commuting Hamiltonians we introduce the Lax matrix L{u) (mon- 
odromy matrix) as a product of n local Lax matrices ii{u) 

L{u)=in{u)...£2{u)ii{u) (2.2) 



(«;Xj,X,)= ( ^^"^^^^ bxi J (2.3) 



where 6, Cj G C are parameters of the model, and u is the so-called spectral parameter 
of the Lax matrix. 

Denoting by id2 the unit 2x2 matrix and introducing notations for the tensor 

1 2 

products i = i ^ id2, i = id2 ®^ one establishes the r-matrix identity |T5[ 

{\i{ui)}j{u2)} = [ri2{ui - U2)}i{ui)ij{u2)]6i„ ri2(n) = -U^u (2.4) 

where 7^12 is the permutation operator in (S> C^. From ( p.4| ) the corresponding 
identity for the monodromy matrix 

{L{ui),Liu2)} = [n2iui-U2),L{u,)L{u2)] (2.5) 



is derived in the standard way ||T5l which, in turn, ensures the commutativity of 
the spectral invariants t{u) and d{u) of the matrix L{u) defined as coefficients of its 
characteristic polynomial 

det{v - L{u)) = v'^ -t{u)v + d{u). (2.6) 

Since det ii{u) = h{u — Cj) the determinant d{u) = detL(n) = nr=i ^("^ ~ Cj) is 
scalar, and the only nontrivial spectral invariant is the trace t{u) 

t{u)=tiL{u) = L^i{u) + L22{u) (2.7) 

which serves as a generating function of commuting independent Hamiltonians Hi 

n 

t{u) = u" + Y,{-iyH,u''-\ (2.8) 

i=l 



As a corollary of (|2.5|) we have the commutativity of t{u) 



{t{m),t{u2)} = (2.9) 

and, consequently, the commutativity {Hi,Hj} = of the Hamiltonians Hi. 
A direct calculation shows that 

n -1 1 " 

H, = N + J2c,, Hl-H--Y^cl (2.10) 

i=i ^ ^ i=\ 

where 

n " /l \ 

N = Y.XiX,, ii = Y.\ + Cia^.X, + fox^+iX^ (2.11) 

j=i j=i ^ 

ensuring that the polynomial ring of commuting Hamiltonians contains the number 
of particles N and the Hamiltonian H. 

Note that the r-matrix ryiiu) in (p.4|) is the same as for the isotropic Heisenberg 
magnetic chain and Toda lattice which puts these integrable models into the same 
class. Indeed, the Toda lattice is a degenerate case of the DST chain. To demonstrate 
it, it is sufficient to make a constant shift u ^ u ^ of the spectral parameter in 
ii{u) given by ( p.3|) and take the limit 

6^0, Xj = e^^(r^ + pj) + 0(6), Xj = e~'?^ (2.12) 

contracting the 'oscillator' algebra {xi,Xi,XiXi) into the Euclidian Lie algebra (e^"^'. 
Pi). In the limit ii{u) turns into the elementary ^-matrix for the Toda lattice: 

(the shifts q become irrelevant since they can be absorbed into a simple canonical 
transformation pi ^ Pi — Cj). On the other hand, the DST model, in turn, is a 
degenerate case of the Heisenberg XXX magnet corresponding to the contraction of 
the su{2) Lie algebra into the oscillator algebra. The DST model holds an intermediate 
place between Heisenberg and Toda models. 

In the present paper we take the Hamiltonian point of view on Backlund tran- 
sformation according to which the Backlund transformation Bx is a one-parameter 
family of simplectic maps from the canonical variables {X, x) to the canonical variables 
(Y,y) possessing certain characteristic properties (see [|l^ for a detailed discussion). 
For Hamiltonian integrable systems allowing a description in terms of the r-matrix 



algebra ( p.5|) there has recently been found an algorithmic method for constructing a 
Backlund transformation [|l6l , |l^. The method having been described in detail in the 
cited papers, we present here only the results. 



As in the case of the periodic Toda lattice [12, O], it is convenient to describe the 



canonical transformation Bx in terms of the generating function 

Fx{y I X) =nX + ± [^^^^^ + (A - q) ln|^^) • (2-14) 
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y dFx 1 A-c, 

Aj = -— = — H (2.15a) 

dxi Vi Xi - hyi+i 

r. = -^ = tx._, + i!ii^i*±ix.. (2.15b) 

oyi Vi 
To prove that Bx preserves the Hamiltonians Hi 

H,{X,x)=H,{Y,y) (2.16) 

we proceed in the same manner as in [|l^, Q for the periodic Toda lattice. Introducing 
the matrices 

one then directly verifies the equality 

Mi{u)ii{u; Xi, X,) = U[u- Yi, y,)M,^^{u) (2.18) 
from which it follows that Bx preserves the spectrum of the Lax matrix L{u) 

M^{u, X)L{u; X, x) = L{u; Y, y)M^{u, A) 



which, in turn, ensures the invariance of t{u) and, therefore, of Hi (|2.16|) . 



To formulate the spectrality property [|13[ of the Backlund transformation we intro- 
duce the quantity /i. canonically conjugated, in a sense, to A 

, dFx " (A - c,)byi A (A - c,)by, 

C'A ~{ byi+i - Xi f-J^ byi+i - Xi 

The spectrality of the Backlund transformation means that the pair (A, /i) lies on 
the spectral curve of the Lax matrix 

det(/i - L(A)) = 0. (2.20) 



To prove it, we again follow DTSf- We observe that for m = A the matrix Mi{u) 
degenerates 

and its null-vector uoi can be found explicitly: 

Mi{X)oj, = 0, u,= ( ^l"-' ) . (2.22) 



Noting then the identity 



(A - Ci)byi 

ii{\)uji^i = Ui, (2.23) 

Oyi+l ~ 
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we conclude that 

L{\)iUn = fXCOn, (2.24) 

whence ( |2.20[ ) follows immediately. 

The commutativity Bx^ o Bx^ = Bx^ o -Bai is an immediate consequence of the 
invariance of Hamiltonians and their completeness, see . 

Note that M~^(m) and ii{u) have, as functions of u, essentially the same structure, 
up to a shift of u and a scalar factor. The fact is by no means a coincidence, see [0 
for a detailed explanation. 



3. Dual Lax matrix 

We conclude the study of the classical case with presenting the dual Lax matrix 
and the dual Backlund transformation for the DST model. In |T0| two different Lax 
matrices were found for the integrable DST system, the 2x2 Lax matrix L{u) and also 
the n X n Lax matrix. This bigger Lax matrix did not contain a spectral parameter. 
Here we present an n x n Lax matrix C{v) containing a spectral parameter v which is 
dual to L{u) in the sense that the corresponding spectral curves are equivalent up to 
interchanging the spectral parameters u and v 

(6" - v) det {u - C{v)) = det{v - L{u)). (3.1) 

To produce the dual Lax matrix C{v) we take an eigenvector 9i{u) of L{u) corre- 
sponding to the eigenvalue v (for brevity, we will not mark the dependence on u in 
0) 

L{u)di = v9i (3.2) 

and define by induction 6i as 

6i+i = ii{u)6i, i = l,...,n. (3.3) 



From (3^) it follows that 9n+i = v6i. The function 6i{u), when properly normal- 
ized, is called Baker- Akhiezer function. Denoting the components of the vector 6i as 
(Pi and ipi we write down ( p. 3D explicitly as 

(3.4) 



fi+i \ ^ ( u- Ci- XiXi bxi \ [ 



Then, splitting the components and taking into account the quasiperiodicity con- 
dition 6n+i = v9i we arrive to the following linear equations for yjj and ipi. 



Uipi = v^i+i + {ci + XiXi)i^i - bxiipi, i = 1, . . . ,n - 1 (3.5a) 

U(pn = V(pi + (C„ + XnXn)iPn " bXnIpn, (3.5b) 
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ilJi+i = -Xi(pi + bilJi, i = l,...,n-l 



(3.6a) 



Vtpi = -XniPn + b^n- (3.6b) 

Eliminating ipi we can write down the linear problem for the vector $ with the 
components ipi in the matrix form: 



C{v)^ = u^, $= ... (3.7) 




where the matrix C{v) defined as 

C(v) = {v- 6")-i E",fc=i b^+^-'^x.XkEjk + vEni 

+ J2j>k V'^XjXkEjk + YJj=i CjEjj + I]"=i Ejj+i (3.8) 

is the dual Lax matrix we were looking for. Here Ejk is the nxn matrix with the only 
non-zero entry {Ejk)jk = 1- The proof of the identity ( p.l| ) is an exercise which we 
leave to the reader. For the case 6=1 and v = —1 the dual Lax matrix for the DST 
model was first found in |TU|. For examples of Lax matrices duality in other integrable 
models see 



The Backlund transformation corresponding to the dual Lax operator C{v) is 
given by the same equations (|2.15|) and ( p.l9| ). The important difference, however, 
is that now yU is a free numerical parameter of Backlund transformation whereas A 
becomes a dynamical variable determined from the equation ( p. 191) . The equality 



( p. 19 ) is now reinterpreted as the equation defining the variable A. The generating 



function of is the Legendre transform of Fx{y \ x) with respect to A. 

The properties of the dual Backlund transformation are proved in the same 
manner as those of Bx (see also |TB| for the Toda lattice case). For the proof we need 
a matrix Ai{v) playing for C{v) the same role that Mn{u) played for L{u). 

Let Oi be defined as 0, = Mi.iOi. From (|]l|) it follows that Oi is a Baker- 
Akhiezer function for ii{u;Yi,yi). The first component of the equality 6i = Mi_i6i 
reads (pi = ipi — byiipi. Substituting ipi from the solution to the system (|3^ ) we obtain 
the correspondence $ = Ai{v)^ with the matrix A^(f) defined as 

n n 

Miv) = iv- 6")-i J2 b'^^'-'yjXkEj, + J2 f^~%X,Ejk + E ^n- (3-9) 

j,k=l j>k j=l 

The invariance of the spectrum of C{v) follows from the identity 

M{v)C{v;X,x) = C{v;Y,y)M{v). (3.10) 
The spectrality is expressed as the identity 

det(A-/:(/i)) =0. (3.11) 
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To prove it is sufficient to notice that the matrix Ai{v) degenerates as f = /i 

detA^(/i) = 0, (3.12) 
and the corresponding null-vector Q defined by the recurrence relation 

^ = ^i^i^^, ^ = l,...,n-l (3.13) 

ili Xi - byi+i 

is, by virtue of ( p.lO| ), also an eigenvector of C{fi) corresponding to the eigenvalue A 

jC{l2)n = Xn. (3.14) 

Since the Toda lattice is a degenerate case of the DST model, the nxn Lax matrix 
for the Toda lattice can be obtained, as one could expect, from our C{v) matrix in 
the limit 6 — >■ 0, as in ( |2.12| ). The result is a variant of the standard nxn Lax matrix 
for the periodic Toda lattice |ll9| : 

£{v) = b-^ + C^^{v) + 0{b), (3.15) 



see 



+ EUiPj + + EU + (3.16) 

Similarly, from M.{v) one obtains the corresponding matrix for the Toda lattice. 



The Poisson brackets for both dual Lax matrices C{v) can be expressed in the 
generalized r-matrix form ^ 



12 1 2 

{C{vi),C{v2)} = [ru{vuV2),C{vi)] - [r2l{Vl,V2), C{v2)], 



the 'non-unitary' r-matrix having the form: 

1 



'^12(^^l,^^2) 



Vi - V2 



kj 



(3.17) 



(3.18) 



^2 I] + t^l I] Ejk 

k>j k<j I 

and r2i(fi,f2) = Vr{v2-,V\)V , where V = Y.]k=i^jk ® E^j is the permutation matrix 
in C" (g) C". 

The non- unitary r-matrix ( 3.18|) in case of Toda's Lax matrix can be unitarized 
by a gauge transformation: 



getting for the new Lax matrix ii(f ) the standard unitary A„_i-type r-matrix 



(3.19) 



ti2(t^i,f2) = z — :7 2^ Ejj 



Vi - V2 



1 



Vi - V2 
1^12(^1,^2) = -^21(^1,^^2) 



k>j k<j 



{h{Vi),l{v2)} = [Xl2{Vl,V2),hvi)+hv2)], 



(3.21) 
(3.22) 



see, for instance, ||T9|| . 
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4. Quantization 



In the quantum case the canonical momenta Xi are replaced with the differen- 
tiations di = d/dxi (having no intent to discuss the conjugation properties of the 
observables, we discard the factor ih to simplify the notation). To preserve the com- 
mutativity of the Hamiltonians Hi upon quantization one needs to choose the operator 
ordering in a special way. 

The necessary algebraic framework is given by the Quantum Inverse Scattering or 
the i?-matrix O, pif] method. Defining the local quantum Lax matrix as 



one verifies the commutation relation 
where 

Rl2{u)=U + Vl2 (4.3) 

is the standard S'Z/(2)-invariant solution to the quantum Yang-Baxter equation. The 
quantum Lax operator, or monodromy matrix, L{u) and its trace t{u) are defined 
then by the same formulas ( p.2|) and (|2.7] ) as in the classical case. From ([4.2|) one then 
derives in a standard way the similar commutation relation 

12 2 1 

Ri2{ui - U2)L{ui)L{u2) = L{u2)L{ui)Ri2{ui - U2) (4.4) 

for L{u) from which the commutativity of t{u) 

[t{m),t{u2)] = Q (4.5) 

follows immediately. The commutative quantum Hamiltonians Hi are defined then, 
like in the classical case ( |2.8| ), as coefficients of the polynomial t{u). It is easy to 
see that Hi is a differential operator of order i leaving invariant the space C[af] of 
polynomials of xi, . . In particular. Hi and H2 are given by the formulas ( |2.10| ) 
with N and H given by (|1.2|) and (|L1|) , respectively. 

The main problem in the quantum case is the spectral problem for commuting 
differential operators, quantum Hamiltonians {ifj}"^^: 

Hii){xi,...,Xn) = hii){Xi,...,Xn), il{xi, . . . , Xn) ^ €.[x\. (4.6) 

One can describe the spectrum and eigenvectors of Hi, or, equivalently, t{u) using 
the well-developed machinery of algebraic Bethe Ansatz [^. Defining the vacuum 
state |0) as the unit function |0) (x) = 1 in C[x] we note that 

L2i|0)=0, Luiu)\0) = aniu)\0), L22(m) |0) = a22(M) |0) , (4.7) 
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where 



a22{u) = 6". 



(4.8) 



i=l 



Defining the Bethe vector ipi;{xi, . . . , Xn) G C[x] parametrized by m complex num- 



bers Vj as 



^PiliXi, ...,Xn) = \Vi,...,Um)= L^ivi) . . . LuiVm) |0) (4.9) 

one can prove [pl[], using the commutation relations (|4.4|), that \vi, . . . ,Vm) is an 



eigenvector of t(u), for any u G C, if and only if the parameters vj satisfy the system 
of algebraic Bethe equations 



and the corresponding eigenvalue r(-u) of t{u) 



t{u) \vi,...,Vm) = r{u) \Vi,..., Vm) 



is given by the formula 



t{u) = an{u} [[ ha22{u}[[ 



j=l U-Vj .^^ U-Vj 



(4.10) 



(4.11) 



(4.12) 



It is usually assumed that Bethe eigenvectors are complete, at least for generic 
values of parameters. The proof of the conjecture is, however, a difficult task, and is 



available only for a few models, see [21] for a discussion 



In his seminal study |11| of the integrable XYZ and XXZ spin chains R J Baxter has 
pointed out that the equations similar to ours ([4.10| ) and ([4.12|) can be reformulated 
equivalently as a finite- difference equation in a certain class of holomorphic functions. 
Adapting his reasoning to our case we introduce the polynomial 0(A; v) in A whose 
zeros are the Bethe parameters Vj: 



0(A;t;) = n(A 



A G C. 



(4.13) 



It is easy to see then that the following finite-difference equation of second order 
for 0(A; v) 

0(A; v)t{X) = aii(A)0(A -l;v) + a22(A)0(A + 1; v) (4.14) 



is equivalent to the system of equations ( |4.10|) for {vj}'jLi and to the equation ([4.12|) 
for r(A). To show this, it is sufficient to divide both sides of ( [1.14| ) by 0(A) and take 



residues at A = Vj. The equation ([4.14|) is called Baxter's or separation equation. 
The reason for the latter name is that an identical equation arises when solving the 



model via the separation of variables method (see for more on relation between 
Q-operator and quantum separation of variables). 

Now we are able to describe the problem we are going to study in the remaining 
sections of this paper. We are looking for a one-parameter family of operators Q\ 
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acting in C[x] such that Q\ shares with t{u) the same set of Bethe eigenvectors, and 
the eigenvalues q{\) of Qx 

Qx \vi, ...,Vm) = g(A) \vi,..., Vm) (4.15) 

are polynomials in A satisfying Baxter's equation ([4.14| ). Up to a normalization coef- 
ficient K^, depending on the eigenvector, the polynomials g(A) are proportional to the 
polynomials 0(A; v) defined by ( [4.131 ): 

g(A) = «:,^0(A; Vi, . . . , Vm) = k^X"^ + 0(A'"-i), A -> oo. (4.16) 

Instead of dealing with eigenvectors and eigenvalues it is more convenient to cha- 
racterize Qx by the following operator identities which are equivalent to the above 
characterization, assuming the completeness of Bethe eigenvectors. We demand that 
Qx commute with t{u) 

[t{u),Qx]=0 (4.17a) 

and self-commute 

[Qa,,Qa,]=0, (4.17b) 
as well as satisfy the finite-difference equation 

n 

Qxt{\) = Qx^i - c.) + &"Qa+i. (4.17c) 

■t=i 

In addition, the eigenvalues of Qx should be polynomial in A 

g(A) G C[A]. (4.17d) 

The above conditions by no means define Qx uniquely. Apparently, one can con- 
struct infinitely many Q-operators just by fixing arbitrary normalization coefficients 
Kij for each eigenvector l^) in (|4.16| ). The difficult problem is to find an explicit ex- 



pression for a Q-operator. Baxter succeeded in solving the problem in case of XYZ 
and XXZ spin chains, having given an expression for Qx as a trace of a monodromy 



matrix However, his formulas do not survive passing to the limiting case of the 
XXX spin chain, governed by the SL{2) invariant R matrix (|4.3| ). 

In the case of quantum periodic Toda lattice, which is another model governed 



by the R matrix ( [4.3| ), a solution was found by Pasquier and Gaudin [jT2|. Instead 
of trying to construct Qx as trace of a monodromy matrix, they considered Qx as an 
integral operator 

Qx : ij{x) ^-^ J dxi . . . J dxn Qxiy \ x) i/j{x) (4.18) 

having given an explicit expression for its kernel Qx{y \ x). They also discovered an 
important relation between the kernel Qx{y \ x) and the generating function Fx{y \ x) 
of the classical Backlund transformation expressed by the semiclassical formula 

Qx{y I x) ~ exp [-^F'xiy | f)) , ^ ^ 0. (4.19) 
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The classical Backlund transformation Bx is thus the classical limit of the similarity 
transformation O ^— QxOQ^^. 



Recently, it was found [|14| how the original Baxter's construction ||TT| can be gener- 
alized to produce Q-operators for the models governed by the Ai-type R matrices, such 
as the XXZ spin chain and sine-Gordon model. According to [jl4|, Qx is constructed 
as the trace of a monodromy matrix built of the local Lax operators corresponding, in 
the auxiliary space, to the special infinite-dimensional representations of the quantum 
group Ug[sl2] (g-osd//ator representations). 

In the subsequent sections we construct a Q-operator for the quantum DST model 
and prove its characteristic properties. Our approach combines those of [12| and of 
Similarly to 



14 



14 



we construct our Q-operator as the trace of a monodromy 
matrix with an infinite-dimensional auxiliary space. In the spirit of [|l^, we find Qx 
as an integral operator acting in C[x] and present several equivalent expressions for 
it. 

The Q-operator being found as an integral operator will give integral equations 
for the eigenfunctions ip^. The advantage of this transformation of the differential 
spectral problem into integral spectral problem is that it gives an alternative to Bethe 
representation of multivariable special functions. The general approach of constructing 
a Q-operator for a given integrable system will be of even greater importance in 
situations when Bethe ansatz does not work. 



5. Construction of the Q-operator 

The structure of Qx is similar to that of t{u) given by (|2.2|) and (|2.7|) . We construct 
Qx as the trace of a monodromy matrix built of the elementary blocks ^x-ci ■ Suppose 
that Ma is a linear operator from C[s,x] into C[t,y]. The spaces C[x] and C[y] are 
referred to as quantum spaces and C[s] and C[t], respectively, as auxiliary ones (see 



np). To construct Qx we introduce n copies ^x-c- '^^ assuming that M^*!^. 



C[si,Xi] C[s,j+i,?/j] (remember the periodicity convention n + 1 = 1) and extending 
^x-c- ^[^j] U 7^ ^) ^ the unit operator. The monodromy matrix 
acts then from C[si, x] into C[si, y\, and Qx is obtained by taking trace in the auxiliary 
space C[si]: 

Qx = tTsMxL.---^^xU- (5.1) 
Supposing Ma to be an integral operator 

Ma : tjj{s,x) J dx J dsTlx{t,y \ s, x) ^{s, x) (5.2) 

we have for the kernel Qx{y \ x) of Qx 

Qxivl X) = j dSn... j dSi [|7^A-c,(s^+l,^/^ | Si.Xi). (5.3) 

i=l 
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To ensure the commutativity \t{u),Qx\ = it is sufficient to demand that Ma 
intertwines 

M{u-X; t, dt)i{u; y, dy)Rx = Rxi{u; x, d,)M{u - A; s, 9,) (5.4) 

the local Lax operator i{u) and some other representation Ai{u — X) of the same 
algebra ( [4.2| ) 

12 2 1 

R{ui - U2)M {ui)M {U2) = M {u2)M {ui)R{ui - U2) (5.5) 

with the same R matrix ( |4.3| ). The proof of (|4.17| a) follows then by a standard argu- 
ment P, p. 

Similarly, to prove [Qai,Qa2] = ( ^4.17| b) it is sufficient to establish the Yang- 
Baxter identity 

J dti ^ dt2 j dy7lx^-\.,{wi,W2 \ ti,t2)Tlx^iti, z \ si,y)Tix^{t2,y \ S2,x) 
= Jdti J dt2fdy 11x2(102, z I t2,y)Tlx^{wi,y \ ti, a;)7^Al-A2(^l, ^2 I ^1,52) (5.6) 

with some kernel Tlx- 

The representation J\4{u — \) of the algebra ( |5.5| ) should be chosen in such a way 
that the resulting Qx, as a function of A, satisfy Baxter's finite-difference equation 
( |4.17| c) and have polynomial eigenvalues ( [4.17| d). As we shall show, for this purpose 
one can take 

M{u;s,d,)=hzf' l) (5-7) 



coinciding essentially with £{u) with 6 = 1 and q = 0. The representation ( |5.7|) plays 
for the Yangian 3^[s/2] the same role as the g-oscillator representation plays for the 
quantum group Uq [s/2] in |14l . Having fixed M. (u) by ( |5.71 ) we get from ( p.4|) a set of 
differential equations for the kernel 7lx(t,y \ s,x) of Ma 

f u + 1 + xd:r bx \ [ u - X + 1 + sds s \ , . . , . 

= [ b )[ 1 j^A(t,yU,x) (5.8) 

(in the right-hand-side we have used integration by parts and the identities d* = —dx, 
{xdxY = ~dxX = — 1 — xdx). The equations ( |5.8|) determine Tlx up to a scalar factor 
Pa: 

nx{t,y \ s,x) = px6{s -by)y~'^exp\^^—^j . (5.9) 



It remains to choose the factor px in ( |5.9| ) and the integration contour in ( |5.2| ) in 
such a way that Ma : C[s, x] C[t, y, A]. 

We describe first the final formula for Ma and equivalent expressions and then prove 
the polynomiality property. As the basic definition of Ma we choose the following 
formula: 

Ma : V(s, x) ^ ^r(A + 1) ^ d^e'^i-^-^-' i^{by, y^ + t). (5.10) 

14 



The infinite integration contour 7 is shown on Figure 0. The branch of the many- 
valued function (— (|5.10|) is fixed by making a cut along (0, 00) and assuming 



that — vr < arg(— ^) < vr. 

7 



Figure 1: Integration contour 7 

From ( |5.10|) it is apparent that M function of A, is analytic in C except the 

poles A = — 1, —2, ... of the factor r(A + 1). As shown below, in fact M.\ continues 
analytically on the whole complex plane. 

Indeed, for Re A < one can pull the contour 7 over the cut (0, 00) and replace 
I-y /(O with /q°° d^ [f{C, — iO) — + iO)] which results in the formula: 

1 r°° 

R^:^^{s,x)^.^^J^ d^e-^r^-'^^iby,y^ + t), ReA<0 (5.11) 



which is analytic in A = —1, —2, .... The branch of ^ ^ in (|5.11 ) is fixed by the 
condition arg ^ = 0. 



To put Ma in the form ( ^.2|) convenient for checking the intertwining relation ( ^.8|) 
one has to make the change of variables x = yC,+t in ( p.lO[ ) . The result is given by the 
formula (|5.2|) with the kernel Tlx given by the expression ( ^.9|) with the scalar factor 
Pa = ^r(A + 1) and integration in x taken over the contour 7' = 1/7 + t. As for the 
integration contour in s, it needs only to pass through the point s = by because of the 
factor S{s — by) in Tlx- 

In the same way, from ( ^.111) one obtains again the formula ( ^.2|) with the kernel 



]9D with the different scalar factor px = l/r(— A) and integration in x taken over 
the ray starting from x = t and going in the direction of y/ \y\. 

Now we have the full description of the operator Ma and can start to study its 
properties. By construction. Ma satisfies the relation ( |5.4| ) from which the commuta- 
tivity [t{u),Qx\ = ( [4.17| a) follows. By direct calculation one can establish also the 



Yang-Baxter identity ( |5.6| ) with the kernel Tlx = Tlx\f,^i proving thus the commuta- 
tivity [Qxi,Q\2] = ( [l.lTj b). The proof of the remaining properties of Qx from the 
list presented in section is given in sections and . 

We conclude this section with giving an alternative description of Ma in terms of 
the polynomial bases which complements the above ones in terms of integral operators. 

To calculate explicitely the action of Ma on the monomial basis s'^x^ in C[s,x] 
one puts ip{s,x) = s^x^ in (|5.10| ), then expands the binomial {y^ + ty and applies, 
termwise, Hankel's integral formula 

/^d5e-^(-0- = -^. (5.12) 
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Using the Pochhammer symbol (c)m = r(c + m)/r(c) = c(c + 1) . . 
one can write down the resuh as 



1 

771=0 



EC) {-^)mt'-"'y"''-'b' = fb%{X;t/y) 



(c + m — 1] 



(5.13) 



where Cm{\b) are the so-called Charlier polynomials 



Cm(A; b) — 2-fo 



The formula ( p.l3| ) proves the polynomiality '■ C[s,x] ^ <C[t,y,X\. Note that 
the normalization of Ma is chosen in such a way that • 1 1- 

The action of Ma on polynomials can be described in an even more compact way. 
Substituting = s'^lx — t)^ into ( p. 10 ) and using again Hankel's formula (|5.12|) 

one obtains the most economic description of Ma 



ix : s\x-ty ^y^+\-X)jb''. 



(5.14) 



In the end of the next section we will discuss the formula (|5.14|) and similar ones 
in more detail. 



6. Analytical properties of the Q-operator 



To produce a description of Qx as an integral operator ( |4.18| ) we substitute the 
expression ( ^.9|) for the kernel Tlx found in the previous section into the formula ( ^.3|) . 
The integration in Si is easily performed due to the delta-function factors in Tlx and, 
corresponding to the two choices of the factor px in ( |5.9|) and integration contour in 
we obtain two equivalent descriptions of Qx- 
The first formula for Qx is given by ( [4.1^ ) with the kernel Qxiy \ x) 



Qxiy\x) = Y[ U7,(A; yi+i, y^, Xi) 



i=l 



where 



Wi{X; yi+i, yi, Xi) = ^T{X + 1 - q) ( — — 



-A-l 



exp 



' byi+i - Xj 
Vi 



(6.1) 



(6.2) 



and integration in Xi is taken over the contour 7^ = yf^ + whereas the contour 7 
is defined in the previous section. 

The alternative formula is given again by ( [4.18| ) with the kernel Q,x{y \ x) 



(6.3) 
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where 

~ / , X yi^ (xi- byi+iV'^^"^ ( bvi+i -Xi\ 

Wi(A; y,+i, yi, Xi) = — — exp (6.4) 

r(Q -X) \ Vi J \ Vi J 

and integration in Xi is taken over the straight ray starting from Xj = 6|/t+i and 
extending to infinity in the direction yi/ \yi\. 

Note that the kernels Qxijj \ x) and QxijJ \ x) satisfy the semiclassical condition 
( [4.19|) which, taking into account our quantization convention — i^ = 1, takes the 



following form (up to insignificant A-depending factors): 

Q\{y I x) ~ exp(-FA(y I x)) (6.5) 



with the generating function of Backlund transformation given by ( |2.14| ). Actually, 
the semiclassical approximation is almost exact, up to a minor quantum correction 
Cj — A I— > c,; — A — 1. This fact supports our thesis on the intermediate place of the DST 
model, with regard to complexity, between the Toda lattice and the generic XXX spin 
chain. For comparison, in the case of the Toda lattice the semiclassical formula for 



Q\{y I x) is plainly exact ||T2[, whereas for the XXX spin chain there is little hope for 
such a simple result. 

For the purposes of the present section we shall need the expressions for Q\ similar 
to the formulas (|5.1CI|) and (|5.11| ) for M^. The corresponding formulas are produced. 



respectively, from ( |6.1| ) and ( |6.3| ) by the change of variables Xi = yi^i + hyi+i. 
The analog of (p.lO|) is the formula 



Qa : i^{x) ^ fd^,... fdu miO • • , yS + by^+u • • •), (6-6) 

where 

" i 

>^a(0 = n ^r(A + 1 - c,)e-^^{-^,y^-'-\ (6.7) 

i=l ^ 

valid for any complex A, except the poles \ = Ci — k, {i = 1, . . . , n; k = 1,2,...) of 
r(A + 1 — Ci). The branch of each of many-valued functions {—^iY^~^~^ in (|6.6| ) is 
fixed by making a cut along (0, oo) and assuming that — vr < arg(— ^j) < tt. 
The analog of (^.11|) is the formula 



QX : ^{X) ^ dCl . . . / dCn >Va(0 ^(. . . , y^^^ + %i+l, • • •), (6-8) 

JO Jo 

with the kernel Wa 

valid for ReA < minReCj. Together, the formulas ( |6.6| ) and ( |6.8| ) define Q\ as a 
holomorphic function of A G C 

In the rest of this section we will show, that Qx maps polynomials in x into 
polynomials in y and A and derive explicit formulas for its action on the monomial 
basis in C[x\. 
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Before considering the general case we will give a brief account of the simplest 
n = 1 case. In this case we have only one variable x = xi, the Lax matrix simplifies 
to L{u) = i{u), so, without loss of generahty, one can put ci = 0. Trace of L{u) gives 
rise to only one integral of motion (number of particles A^) 

t{u) =tTL{u) =u- N + b, N = xd. (6.10) 

We assume that N acts in the space C[x] of polynomials of x spanned by the 
eigenbasis {x"^}'^^q of 

N : x"" ^ mx"", m = 0,1,2, ... . (6.11) 
For n = 1 and c = the formula ( |6.6| ) defining the Q operator turns into 

Qx : ^(x) ^ ^r(A + 1) ^ d^e-^i-^-'-' ^{y{^ + b)), A ^ -1, -2, . . . (6.12) 

and ( |6.8| ), respectively, into 

1 

^ ^(^) ^ fPAy i, dee-«r'"V(2/(e + &)), ReA<0. (6.13) 



Similarly, from ( |6.1|) and ( |6.3|) one gets, respectively, 

y, 



Q^:^{x)^'^T{X + l) f dxy-'lb--] e-^/y^{x), 7' = 1/(7 + 6) (6.14) 
2n Jy \ y J 

and 



e r°° , 1 / X 



Qa:^(x)^— --/ dxy-U--b] e-^/^^(x), y > 0. (6.15) 
r(-A) Jby \y J 

To calculate explicitly the action of Qx on the basis x"^ one puts iIj{x) = in 
( |6.12|) , then expands the binomial (^ + b)"^ and applies, termwise, Hankel's integral 
formula (|5.12| ). This calculation is very similar to the calculation of M.xs'^x^ given by 
the formula (|5.13| ). The result is that the monomials {x™}5^^o the eigenvectors of 
Qx: 

Qa : X™ ^ gm(A)|/™, (6.16) 

the corresponding eigenvalues being polynomials in A of degree m expressed in 

terms of the Charlier polynomials Cm(A; b) as 

QmW = E ■ (-A),6™-^ = 6"^C„(A; b) (6.17) 



j=0 



J 



(compare to ( ^.13| )). 

As an immediate consequence, we have the commutativity [Qx, N] = 0, as well as 
( [4.17| a) and ([4.17| b). Another corollary is that Qx maps C[x] into C[y, A]. Note that 
formula ( 5.17 ) implies the normalization Qx : 1 ^ 1. 
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One can use the integral operator Qx to derive few well known formulas for the 
orthogonal Charlier polynomials. For instance, putting iIj{x) = x™" and ?/ = 1 in (|6.12|) 
or in ( |6.13|) one obtains integral representations for Charlier polynomials: 



C^(A;6) = ^r(A + l) J d^e^i-^)-^-' (l + 



(6.18) 



and, respectively. 



r(-A) Jo 



d^e-^r'^Ml 



b 



(6.19) 



see 



H). 



Equating '?/'(x) in (|6.12| ) or (|6.13|) with the generating function e*^ = J2m=o x'^t"^ / t^^- 
of the monomials and taking the integral one gets the generating function of 
Charlier polynomials 



oo ±m 



m=0 



ml 



(6.20) 



The recurrence relation for the Charlier polynomials |2^ is equivalent to the 
finite-difference equation for the polynomials qi{X) 

{X-1 + b)qi{\) = hqi{\ + 1) + \qi{\ - 1), (6.21) 

which coincides with Baxter's equation ( |4.14D for n = 1 and proves, for = 1, the 
operator relation ([4.17| c). 

From the explicit expression (|6.171 ) for the polynomials qm{X) we conclude that 
they are normalized by the condition gm(0) = 6"*, or, alternatively, gm(A) = (— A)*" + 
0(A'"~^), as A — >^ oo. In terms of the operator Qx it is equivalent to 



(see (16.101) for the definition of A^) and, respectively, to 

gA = (-A)^ + 0(A^-i). 



(6.22) 



(6.23) 



The generalization of the above results to the multivariable case is quite straight- 



mi 



forward. To calculate explicitly the action of Qx on the monomial basis x" 
in C[x] one substitutes ipi^x) = . . . x^" into ( |6.6|) , then expands the binomials 
iVi^i + hli+i)^' and uses termwise Hankel's integral formula ( |5.12| ). Recalling the 
definition ( |6.17|) of Charlier polynomials, one obtains the following expression: 



Qx-x\ 



(6.24) 



i=l 



from which it follows immediately that the normalization condition Qx : 1 ^ 1 holds 
and that Qx maps C[x] into C[^, A]. The polynomiality of matrix elements of Qx 
combined with the commutativity [QxiyQxi] (|4.17| b) proves the polynomiality ( [4.17] d) 
of the eigenvalues of Qx- 
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The formula ( |6.24D also allows to determine the normalization ( f4.16| ) of the eigen- 



values of Qx. Taking the limit A — oo in (|6.24|) and using the asymptotics Cm(A; h) 



\/b)"^ + 0(A'" ^) we conclude that, as in the n = 1 case, Qx has the asymptotics 
^231) with the operator given by (|1.2|). In contrast, the equality ( |6.22| ), generally 



speaking, cannot be generalized to n > 1, with the exception of the homogeneous 
chain case Ci = 0, i = 1, . . . ,n, when it is replaced by 

Qo = b^U (6.25) 

where U is the translation operator U : Xi Vi+i- 

As a final remark of this section, we point out yet another way of writing down 
the action of Qx- Substituting ip{x) in (|6.6|) with the polynomials G C[x\ 



i=l 



parametrized by the multi-index m = (mi, . . . , m„) and a vector y = (yi, . . . , yn) we 
obtain, after performing the integrations, an elegant formula for the action of Qx on 



, ,m. 



Qx : n(^^ - by^+T' ^ n(c. - >^)m.yr- (6.26) 

i=l 1=1 

The formula (|6.26|) seems to provide the most compact way to encode the action 
of Qx on polynomials (compare to the formula ( |5.14| ) for the action of Ma)- Some 



caution is necessary, however, when using it, since the parameters y in o;^ coincide 
with the variables in the target space of Qx- One way of interpreting ( |6.26| ) is to 
consider its left-hand-side as a shorthand notation for [Qxo~i'^]z=y- Another possibility 
is to extend the operator Qx onto the polynomial ring C[x, assuming that it acts 
on y trivially: Qx{;iIj{x)^{ii)) = (p{y)Qx{'ip{x))- Formulas similar to ( |6.26| ) also arise in 



the separation of variables for Macdonald polynomials ||2^ . 

It is a challenging problem to take the formulas ( p.l4| ) and ( |6.26| ) as definitions of 
Ma and Qx, respectively, and to build the theory of Qx in an entirely algebraic way. 



7. Baxter's equation 

In the previous sections we have proved all the properties of Qx from the list given 
in section except Baxter's difference equation ( |4.17| c). In this section we give a proof 
of the identity ( [4.17| 3) based on the ideas of |ll2 |. 



The best suited for our purposes realization of Qx is that given by the formulas 
( [4.18|) and ( |6.1| ). Recalling that t{u) = ti L{u) and that L[u) is a 2 x 2 matrix whose 



entries are differential operators in Xj, we can transform the left-hand-side of ( [4.17| c) 
as follows 

[Qxt{Xmy) = tT[QxL{X)ij]{y) =tT [ dx^ Qx{y \ x)L{\)ij{x). 
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Performing integration by parts we obtain 

[Qxt{X)^|J]{y) = tr / dx" [L*iX)Q,iy \ (7.1^ 



where L*{X) is the matrix composed of adjoint differential operators [Lj^)* = L*f^. For 
example, d* = —d, (xd)* = —dx = —xd — 1, and so on. 

Using the factorization ( p.2| ) of L{X) into the product of elementary Lax matrices 
£j(A) and the factorization ( |6.1| ) of the kernel Qx{y \ x) into the factors Wi ( |6^) we 
can represent the kernel of the integral operator Qxt{X) as 

n 

[Qxt{X)]{y\x) = trC(A) . ..il{X)Y[w, = tr(C(A)^„) . . . (^t(A)^i), (7.2) 

1=1 

where 

i*{X) = 1^ ^ ~ + ^ . (7.3) 

The possibility of the factorization (|7^ ) of [QxL{X)]{y \ x) depends crucially on 
the fact that the factors Wi ( |6.2| ) depend each only on one variable Xi. That is why 
we take the left-hand-side of (|4.17| c) to be Qxt{X) rather than t{X)Qx- 

The rest of the proof parallels that of the spectrality property for the classical case 
given in section . Introducing matrices ii{X) by the equality i*{X)wi = Wiii{X) and 
noting that 

d^^ lnwi{yi+i, yi, Xi) = . (7.4) 

Xi - byi+i yi 

we obtain 



X , 1 , 1 U \ I - A - l)y,+l Xi r 

, / A — Ci + 1 + XiOx inWi bXi \ ( r—hi-,^ Vi * 1 /^, ^\ 

Xi - byi+i Vi 



and 

[QxtiX)]iy\x) = Qxiy\x) tr 4(A) . . . £i(A) = Qx{y\x) trZ(A). (7.6) 

We note then that the gauge transformation £j(A) t-^ A'jT,_\£j(A)A^j with the gauge 
matrix 

A'. = ( ; "f ) p'^) 

leaves trL(A) invariant while making ii{X) and, consequently, ^^(A) triangular: 

Nr^\UX)N, = 



(7.8) 




w,{X) J 



where we used the identities 

Wi{X + l) _ - A - l)yi Wi{X -1) _ Xj- byj+i 

Wi{X) Xi-byi+i ' Wi{X) {ci-X)yi' 

As a result, we get the equahty 



(7.9) 



which, obviously, proves ( [4.17| c). 



8. Discussion 



On the example of the quantum integrable DST model we have shown that the 



construction of Q-operator as an integral operator, in the style of the paper |]T2[, and as 
the trace of a monodromy matrix with a special representation of the quantum group 
corresponding to the auxiliary space, in the style of papers |]T3, can be combined 



naturally within an unified approach. The same approach can be applied to other 
integrable models which are more general than DST model, such as the generic XXX 
magnetic chain. This work is in progress and the results will be reported in a separate 
paper. For a particular case of the homogeneous XXX chain a Q-operator was recently 
constructed in p5[] . 

Another interesting problem is to build the theory of Q-operator in a purely alge- 
braic manner starting from the formulas ( |5.14| ) and ( [6.261 ). 



In it is argued that for the models governed by the Ai type R matrices there ex- 
ist two different Q-operators corresponding to two different g-oscilator representations 
of Wg[s/2]. Their eigenvalues correspond, respectively, to two linearly independent so- 



lutions of Baxter's difference equations analogous to ( [4.141 ). In the case of DST model 



the second Q-operator can be obtained if we choose in the formula ( [5.4[ ) another 
representation M.{u — X) of the algebra (|5.5[ ), namely A^(m — A) ~ —A4~^{X — u) 



MW;.,a.)={:g^ J^g]. (8,i; 



The corresponding Q-opsrator has, however, more complex nature than the one 
studied in the present paper. Its eigenvalues, for example, are not polynomial in A. 
The problem is currently under study. 

We can point out the following application of our results to the theory of special 
functions of many variables. Notice that the eigenfunctions of the quantum DST 
Hamiltonians are multivariable polynomials. The obtained family of integral equations 
for those polynomials provided by the QA-operator supplements their representation 
as Bethe vectors and can be used in efficient numerical calculations of these special 
functions, for instance, solving integral equations by iterations. Simple considerations 
of the n = 2 case show that we deal with the discrete multivariable analogues of 
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the Heun polynomials. For special functions of such complexity the found integral 
equations might be the only explicit representations to exist because there is no hope 
to get, for instance, an integral representation. So, the found integral equations for 
the special functions which were initially defined as eigenfunctions of the commuting 
differential operators can be used first, as already remarked, for generating advanced 
numerical methods of their calculation, and, secondly, for finding various asymptotics. 
These applications are being worked on. 
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